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We derive a hypervirial relation for the positronium ion, a three-body bound state of two electrons 
and a positron. It relates expectation values of three operators and reconciles two recently published, 
seemingly inequivalent formulas for the magnetic moment of the ion. As a by-product, the precision 
of the leading binding correction is improved. 
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The positronium ion is no longer an exotic system. This very weakly bound state of two electrons and a positron 
can now be copiously produced with a recently discovered method (TJ, [ 2 | ; up to two percent of positrons can be 
converted into ions Ps _ . This advance stimulates theoretical studies of the properties of the ion, including its magnetic 
moment. The two electrons in Ps - form a spin singlet so the entire magnetic interaction is due to the positron, whose 
gyromagnetic ratio is however slightly modified by the binding. This binding effect, of the order of the fine structure 
constant squared, a 2 , has recently been determined using two different approaches a a , the latter one more precise. 
The result is the expectation value of a combination of three operators, involving the kinetic and the potential energy 
of the positron, and a correlation of position vectors of the three particles. 

The two studies agree on the total numerical value, but find quite different relative sizes of the three contributions. 
It has been thought that the two expressions are not equivalent [J], casting doubt on the correctness of at least one of 
the approaches. Here we derive an identity that relates the three expectation values and thus reconciles both results. 

The identity is of the hypervirial type [5j, [1] ■ In addition to revealing the equivalence of the two published results 
sa, it allows us to eliminate one of the three expectation values and thus present a new expression, simpler and 
more precise, for the magnetic moment of Ps _ . 

A challenge in studying Ps _ is that the wave function of this three-body system is not known analytically. It is 
computed with the variational method, using the nonrelativistic approximation. In order to present the Hamiltonian, 
we denote the coordinates of the two electrons by fj and r* 2 , and use r *3 for the positron. The internal dynamics 
depends only on the relative coordinates r \2 = r? — rj and ri 3 = r *3 — rj, and not on the position of the center of mass 
R = | (fq + r *2 + t%). The momenta of the three particles become 
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and the Hamiltonian, written in atomic units (and with Vf,., = Vij) 
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The bound positron g-factor is the sum of the free particle value <?f, 
three expectation values needed to express Abound are 


1 1 

ri3 7-23 ' 

ee and the binding correction Abound- The 


A = (pi) = -<V? 3 ) = 0.257 532 962 
B = (y-j = 0.339821023 

C = ( j 13 ^/ 23 ^ = 0.046 478 421, 
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where the numerical values, taken from Q, are presented here only to allow the reader to check the relation among 
them, derived below. The values of A and B agree with the more precise evaluations by Frolov [7]. The results for 
Agbound found in the recent studies are (for now neglecting self-interaction corrections, taking gf lee -a 2 ) 


Ref. g 2^1 = - 1 (B - C), 

a z 9 3 

Ref. a ^ = -A-^B+^C. 
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In order to show the equivalence of these expressions, consider the expectation value of the commutator of the 
Hamiltonian (0) with the operator r *23 • Vi 3 . Taken in a stationary state, such an expectation value must vanish, since 
it expresses the change with time of a time-independent operator [ 8 !j. Evaluating the commutator we find the identity 
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In order to eliminate the scalar product of two momenta, we note that both electrons have equal average momentum; 
using m and (0 we find 


(v 1 2 -v 13 ) = ~(v? 3 ), 


( 11 ) 


in agreement with |7], in whose Table II the values of (pi ■ p 2 ) and (pi ■ p^) should have minus signs |9]. Also equal 
are the average potential energy of each electron interacting with the positron, 


— ) = ( — 

J 23 / \ r 13 , 

Substituting the last two equalities into (fl()l) . we obtain the main result 
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It turns out that the two published expressions (0) and © differ by 4/27 times this combination, that is by zero. An 
evaluation of C using (THU) and the precise values of A and B |7] gives a value that differs from (JT]) only in the last 
digit (it is 0 instead of 1). 

Using (fT51) we can eliminate C in terms of A and B. The binding correction becomes 
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where we have included the actual value of the free-positron ^-factor, using the result of Q (rather than approximating 
gfree —> 2 as in i)- Taking the recently measured <?f ree of the electro njl0| and the values of A and B from Q, we 
confirm and improve the result of Q, as well as the less precise result [ij 

Abound = -0.510 551028187 6(6)a 2 + O (a 4 ) . (15) 

The error in the coefficient of a 2 is dominated by the comparison of ^-factors of free electrons and positrons, known 
to be equal to better than three parts per trillion 0- Of course, providing that coefficient with more than four 
or five decimal places is at present only of academic interest, because of the unknown a 4 effects. Nevertheless, the 
improvement of precision shows the power of the hypervirial relation that eliminates a lesser known expectation value 
C in favor of A and B, very well known from the determination of the binding energy of Ps _ . 

The total g-factor of the bound positron is (with a from fl2j or from a comparison of the measured gf re e with QED 

0 ) 


g Ps - = ghee + Abound = 2.002 292117 (3), (16) 

where the error is estimated by ±od Q 1 . 

To summarize, we have demonstrated that the two recently published formulas for the magnetic moment of Ps _ are 
equivalent. Eq. m is but one example of hypervirial relations that can be used to check the accuracy of variational 
calculations. This is a welcome development, now that Ps” has become more accessible to precise measurements 0. 
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